Zero-temperature criticality in the two-dimensional gauge glass model 
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The zero-temperature critical state of the two-dimensional gauge glass model is investigated. It is 
found that low-energy vortex configurations afford a simple description in terms of gapless, weakly 
interacting vortex-antivortex pair excitations. A linear dielectric screening calculation is presented 
in a renormalization group setting that yields a power-law decay of spin-wave stiffness with distance. 
These properties are in agreement with low-temperature specific heat and spin-glass susceptibility 
data obtained in large-scale multi-canonical Monte Carlo simulations. 
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Rigidity is a key concept that unifies description of 
long-ranged order in magnets, superfluids, superconduc- 
tors, and crystalline solids. For two-dimensional (2D) 
systems with an 0(2) symmetry, the helicity modu- 
lus provides a quantitative measure of rigidity, and as- 
sumes a finite value in a pure sample at sufficiently low 
temperatures. 2j 3] There are however a number of exper- 
imental situations, notably granular thin films of high-Tc 
superconductors, where one has to deal with quenched 
disorder that introduces random frustration and ground 
state vortices into the system. 4, 5, 6] Understanding how 
rigidity weakens or disappears altogether is crucial for 
interpreting the often complex equilibrium and dynamic 
behavior of these systems at low temperatures. 7] 

The gauge glass modelj^Q in two dimensions offers a 
good example where the nature of complexity and glassi- 
ness associated with strong disorder can be examined in 
quantitative detail. Two competing scenarios regarding 
the low-temperature phase diagram have been proposed: 
i) The system is disordered at any temperature T > 0, 
but critical at T = 0, as characterized by a power-law di- 
verging glass correlation length fq — T~'^ with a rather 
large exponent i/~2.5. HEIllElESii) There is a 
glass phase of frozen vortices that undergoes a continuous 
transition at Tc ~ 0.2 J. ^3^0,0] Direct measurement 
of rigidity in terms of spin-glass susceptibility has not 
been able to differentiate the two scenarios unambigu- 
ously due to the relatively small system sizes examined 
in simulations. The domain-wall energy analysis'^ has 
provided some insight on the characteristics of low-lying 
excitations, although their true identity remains myste- 
rious. 

In this paper, we describe results of extensive numer- 
ical investigation of the 2D gauge glass model and the 
related Coulomb gas problem. We show that, despite 
strong correlation in vortex positions and the collective 
nature of the ground state, low energy excited states of 
this model take the form of a dilute gas of fermion-like 
vortex-antivortex pairs which have a gapless, excitation 
spectrum. Based on the existence of such elementary ex- 
citations, we present a phenomenological renormalization 



group (RG) calculation that explains the power-law de- 
cay of spin-wave stiffness J{R) ^ with distance R in 
the ground state. Evidence for gapless pair excitations 
is seen in the linear specific heat at low temperatures. 
The stiffness exponent 6 is estimated from the strength 
of the screened disorder potential in the ground state 
of the Coulomb gas problem, and from the spin-glass 
susceptibility data obtained in large-scale multicanoni- 
cal Monte Carlo simulations of the gauge glass model. In 
both cases, we obtain 6 ~ —0.45, in approximate agree- 
ment with earlier studies. 

The gauge glass model is defined by the Hamiltonian, 
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where (pi is the phase variable on site i, J is the coupling 
constant, and Ay are quenched random variables uni- 
formly distributed on [— tt, tt). Summation is over near- 
east neighbor pairs of sites, here on a square lattice. De- 
composing the lattice gradient of the phase field (pt into 
rotation-free and divergence- free components, one arrives 
at the following Coulomb gas model of vortices. [l^ 
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Here is the vortex charge on site i, Ec is the vor- 
tex core energy, a is the lattice constant, and is 
the distance between sites i and j. The random poten- 
tial Vi can be expressed as a sum of dipolar potentials 
produced by the random phase shifts Aij. Its statis- 
tics is specified by a logarithmically growing variance 
(Yi) ~ 27rcrJ^ ln(L/a), and logarithmic spatial correla- 
tions {{Vi — VjY) ~ AiTa.P \w{rijla). For the gauge glass 
model, we have a = 7r^/3 and — 5 J, such that ap- 
proximately one third of the sites are occupied by either 
a vortex or an antivortex in the ground state. 

Due to the high density of vortices and antivortices, 
the RG analysis 18, 19] developed for a dilute Coulomb 
gas can not be applied here. We have investigated 
numerically the minimal energy states of the Coulomb 
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gas Hamiltonian (0) using a greedy algorithm, restrict- 
ing vortex charges to rrii = 0, ±1. Such states are 
constructed through successive addition and removal of 
vortex-antivortex pairs, each time picking a pair state 
of the lowest energy. (Note that moving a vortex to a 
vacant site can also be realized by such a move.) The 
process terminates when the system energy can not be 
lowered further through single pair addition or removal. 
In an improved ground state search, we allow the pro- 
cess to continue until a previous vortex configuration is 
revisited. Our experience shows that, for systems up to 
a linear size L = 16, the algorithm is able to find the 
ground state typically in 100 or fewer trials of random 
initial conditions. 

Figure 1 illustrates results of our finding for a 16 x 16 
system. The ground state vortex (solid circle) and an- 
tivortex (open circle) configuration {m^} is shown in Fig. 
1(a). As expected, there is no apparent order of the vor- 
tices and antivortices. The same is true also for other 
low energy states obtained under the greedy algorithm. 
However, when the excess charge {rhi = rrii — to^} of 
these states is plotted, as shown in Figs. l(b)-(d), we 
see clearly that the low energy states overlap strongly 
with the ground state. The excess charges can be iden- 
tified as elementary excitations that include single vor- 
tex/antivortex movement to a new position (e.g., vortex 
at site A to site B, vortex at site C to a neighboring 
site, antivortex from site G to a neighboring site, and 
antivortex from site H to site I), and insertion (pair D) 
or deletions (pairs E and F) of vortex-antivortex pairs. 

Repeated numerical experiments of this type led to the 
following general observations: i) The ground state con- 
tains a high density of vortex-antivortex pairs of all sizes, 
ii) Low energy excited states can be described as a dilute 
gas of excess vortices from the ground state, iii) Most 
of the excess vortices in ii) form closely bound vortex- 
antivortex pairs. However, pairs of larger size [e.g., the 
vortex-antivortex pair at H and I in Fig. 2(d)] may also 
be present. They usually form a complex with smaller 
pairs that provide screening of the Coulomb interaction. 

To characterize the energetics of low-energy pair exci- 
tations, we rewrite Eq. (0) in terms of the excess vortex 
charge rhi = rrii ~ n^i, 
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Here Eq is the energy of the ground state {m"}, and Vi is 
an effective potential at site i that includes the original 
disorder potential Vi as well as the potential produced by 
the ground state vortex population {m-*}. 

In the original model, different Fourier components 
of the disorder potential are gaussian distributed and 
statistically independent of each other. Their vari- 
ance is given by (y(k)F(-k)) = Na(27rJ)^G{k), where 
iV = is the total number of sites and G(k) = 
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FIG. 1: (a) Vortices (solid circles) and anti-vortices (open 
circles) in the ground state, (b)-(d) Excited states shown in 
terms of excess vortices (solid circles) and antivortices (open 
circles) from the ground state. Energies are in units of J. 



l/[4sin^(fc3;/2) -|-4sin^(fcy/2)] is the lattice Green's func- 
tion. Figure 2(a) shows {V{k)V{-k)) /G{k) at ky = 
and kx = 27r/L, Att/L, . . . , tt, for ground states obtained 
under the greedy algorithm and averaged over many dis- 
order realizations. All energies are in units of J. Statis- 
tical errors are smaller or comparable to symbol size. It 
is evident that, unlike the "bare" disorder potential Vi, 
long wavelength components of Vi are much reduced in 
strength by the ground state vortices. In fact, the data 
can be fitted well by assigning a k-dependent coupling 
constant J{kx,0) ~ | sin(A:2;/2)|~^, with 6 = -0.45, as 
illustrated by the dashed line in the figure. 

We have also examined the size-dependence of the gap 
A between the lowest and second lowest values of Vi in a 
given disorder realization. This quantity offers a measure 
of the spectral density of Vi close to its minimum value, 
which governs vortex movement on scale L. Figure 2(b) 
shows the distribution of A for three different system 
sizes, which is Poisson like. The inset shows a scaling 
plot of P{A), assuming a power-law width A^, ~ L^'^-'^^. 
A nearly perfect data collapse is seen. 

To explain the aforementioned behavior, we consider 
a dielectric screening model of vortex-antivortex pairs in 
a phenomenological RG setting. Let R be the running 
length scale in an RG scheme applied to (jSJ, and J{R) 
be the renormalized interaction strength after pairs of 
size less than R have been integrated out. In an infinite 
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FIG. 2: (a) Scaled power spectrum of the effective disorder po- 
tential against kx for three different system sizes. The dashed 
line indicates a power-law fit with an exponent 2\9\ — 0.9. (b) 
Distribution of the gap energy between the lowest and second 
lowest values of the effective potential. Inset shows the scaling 
plot of the distribution with a stiffness exponent 6 — —0.45. 

system, the energy e of pair states in the size range R to 
R+dR forms a continuum that, for the gauge glass model, 
extends over both positive and negative values. Although 
filling of states with e < is a collective process, the pair- 
pair Coulomb interaction Vpp ^ r^'^ is not strong enough 
to generate a gap at the "Fermi level" ep = O.'20l| Conse- 
quently, we expect treating pairs as independent ferniions 
is qualitatively correct. An external perturbing field E 
shifts the energy of state i by an amount — E • pi, where 
Pi = — r~ is the dipole moment of the pair in state i. 
In the neighborhood of ep, occupied and unoccupied pair 
states may switch under the perturbation. The induced 
polarization is easily calculated, from which one obtains 
the dielectric susceptibility for such a medium at T = 0, 

Xfl = R^Pr{0) X {dR/R). (4) 

Here pr{0) x (dR/R) is the density of pair states per unit 
area at e = 0. Let dl = dR/R, change in J{R) due to this 
group of vortex-antivortex pairs can be expressed as|l9l|. 

dJ-^/dl = 47r2p(0). (5) 

where we have introduced the rescaled density of states 
p{0) = i?^pfl(0). Since J{R) is the only energy scale at 
zero temperature, we may write, on dimensional grounds, 

m = (6) 

where c is some constant, jl^l Combining Eqs. lO and 
© , we obtain a power-law decay of the coupling constant 
J{R) = JbR^ in the ground state, with = —Aii^c. At a 




FIG. 3: Specific heat versus temperature for a 48 x 48 system. 
Dashed line shows a quadratic fit to the low temperature data. 

finite temperature T, Eq. jnj must be modified to take 
into account thermally excited pairs when one reaches a 
scale £,g{T) set by J(^g) — T. This yields a correlation 
length ^ T~^/I^l beyond which rigidity disappears. 

We have carried out extensive simulations of the gauge- 
glass model using a multi-canonical Monte Carlo sam- 
pling scheme that allows one to equilibrate systems of 
size up to L = 48 down to T = 0.05. Computation were 
performed on a PC cluster using up to 64 nodes. The 
parallel setup is ideal for calculating quantities involving 
the overlap q = exp[i((/)" — of configurations 

from two different replicas a and h under the same disor- 
der. Details of our simulation algorithm will be reported 
elsewhere. 

Figure 3 shows the specific heat data for the L = 48 
system, averaged over 26 disorder realizations. Compar- 
ison with data at smaller L (not shown) indicates that 
no significant finite-size correction is present. Data on 
the low temperature side can be fitted to the formula 
cy = 0.5 + 0.32r + 0.8T2, which has a simple physical 
interpretation. The constant term is due to spin-wave 
contributions. The linear part confirms presence of gap- 
less vortex-antivortex pair excitations in the system. 

We have also measured the k-dependent spin-glass sus- 
ceptibility 

XSG (k) = ^CsG( Tij) exp(zk • Vij 

), (7) 





where CsG{i~ij) — (|e*('^»~'^3) | } is the correlation func- 
tion of the glass order parameter. Here the overline bar 
denotes thermal average and (•) denotes average over the 
disorder. Figure 4(a) shows xsg(O) = A^(|gp) against T 
for eight different system sizes. In agreement with pre- 
vious studies, the glass phase correlation grows rapidly 
at low temperatures. For the largest system shown at 
L = 48, crossover to finite-size dominated regime takes 
place already at T ~ 0.3. 
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FIG. 4: (a) Spin-glass susceptibility versus temperature for 
various system sizes L. (b) Temperature derivative of the 
susceptibility ratio at |k| = and 2n/L versus T/J{L) — 
TL~^ for six sizes ranging from L = 6 to 32. See text for 
details. 



Since spin-wave fluctuation introduces a power- 
law decay of Csg('') at short distances with a 
temperature-dependent exponent il{T\ extracting £,g{T) 
from the finite-size scaling ansatzjl^ xsG(k, T, L) — 
L^~^x{L^,^gI L) is somewhat ambiguous. Instead, one 
may consider the ratio a{T) = xsG(|k| = 27r/i)/xsG(0) 
which tends to 1 for ^g ^ and for ^g 3> L. Crossover 
between the two regimes takes place ai — L. Fig- 
ure 4(b) shows the derivative da{T)/dh\T for six dif- 
ferent system sizes ranging from i = 6 to 32. The 
peak position of these curves yields a quantitative mea- 
sure of ^G- Plotting the data against the scaling vari- 
able T/J{L) = T/L^, we obtain the best collapse at 
= —0.45, in agreement with the zero temperature anal- 
ysis presented above. 

To summarize, the picture that emerges from our study 
of the 2D gauge glass model is a ground state that sup- 
ports gapless, fermion-like vortex-antivortex pair excita- 
tions. The energy scale of these excitations decays as 
a power-law of pair size — a consequence of dielectric 
screening by pairs of smaller size. The hierarchical orga- 
nization of low energy vortices, with diminishing energy 
scales on increasing length scales, offers a new paradigm 
of zero-temperature criticality. This description is fully 
consistent with previous numerical work as well as our 
low-temperature specific heat and spin-glass susceptibil- 
ity data obtained under a multicanonical Monte Carlo 
sampling scheme. 

The composite nature of the low energy vortex- 
antivortex pairs of large size, with smaller and localized 
pairs providing screening of the long-ranged Coulomb po- 
tential, suggests possible glassy dynamic behavior. In a 
nonequilibrium context, the effective energy scales asso- 
ciated with long-distance vortex motion may be much 
higher than the equilibrium ones that require screening 



vortices to find their ideal positions. Previous simu- 
lation work on the relaxational dynamics of the gauge 
glass model indicates a rich behavior with multiple size- 
dependent time scales. 0,^3 I* would be interesting to 
re-analyze the data in light of the new insight on the 
energetics of elementary vortex excitations. 
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